Measurement of b hadron masses in exclusive J/tp decays with the CDF detector 
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We measure the masses of b hadrons in exclusively reconstructed final states containing a J/ij) — > 
fl~ decay using 220 pb _1 of data collected by the CDF II experiment. We find: 
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The measurements of the B° , A° mass, m(B®) — m(B°) and m(A°) — m(B°) mass difference are of 
better precision than the current world averages. 



PACS numbers: 13.25.Hw, 13.30.Eg, 14.40.Nd, 14.20.Mr 



In the Standard Model, Quantum Chromodynamics 
(QCD) is the theory describing strong interactions be- 
tween objects of color charge y|. Hadrons are colorless 
particles made up of strongly-interacting constituents. 



Hadron masses are fundamental physical observables and 
their study is the spectroscopy of quark systems bound 
by QCD. The heaviest known hadrons contain a bot- 
tom or b quark 0|. Lattice QCD calculations predict the 
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hadron mass spectra from first principles. Comparing 
these calculations with experimental data is an essential 
test of QCD. Recent advances in lattice QCD, using un- 
quenched methods, will soon allow mass predictions with 
close to the experimental accuracy Q and we look for- 
ward to these new calculations to compare with our re- 
sults. The most precise predictions are those of b hadron 
mass differences. We present here the most precise indi- 
vidual measurements, to date, for the masses of B + , B°, 
B° and A°. 

The data used in this analysis were obtained with the 
Collider Detector at Fermilab (CDF II) operating at the 
■s/s = 1.96 TeV Tevatron pp-collidcr. The data were 
collected between February 2002 and August 2003 and 
corresponds to an integrated luminosity of 220 pb _1 . 
The CDF II detector is described in detail elsewhere 0. 
This analysis relies on the tracking system and the muon 
detectors. The tracking system is comprised of a sili- 
con micro strip vertex detector (SVX II) 5] and a drift 
chamber operating in a 1.4 T solenoidal magnetic field. 
The SVX II system consists of 5 concentric silicon lay- 
ers made of double-sided silicon covering the radii from 
2.5 cm to 10.6 cm. The impact parameter resolution is 
about 40 /im, including a 30 /im contribution from the 
beamspot. The Central Outer Tracker (COT) Q is an 
open cell drift chamber measuring 310 cm in length, with 
an inner radius of 41 cm extending to a radius of 138 cm 
and provides a large lever arm for curvature measure- 
ments. Each cell contains a plane of 12 sense wires tilted 
by 35° with respect to the radial direction to compen- 
sate for the drift Lorentz-angle. The COT is segmented 
radially into eight superlayers. For superlayers 1, 3, 5 
and 7 wires form a ±2° stereo angle with respect to the 
beam direction, while for superlayers 2, 4, 6 and 8 wires 
are oriented along the beam direction. The measured 
momentum resolution is ct{pt)/pt ~ 0.15%pr/(GeV/c). 
Muon detectors consist of multi-layer drift chambers lo- 
cated around the outside of the calorimeters 0. The 
central muon system covers a range in pseudorapidity of 
< 0.6. The central muon extension extends the pseu- 
dorapidity range to 0.6 < \rj\ < 1.0. 

Data are selected with a three-level trigger system. 
The Level 1 portion of the dimuon trigger uses the ex- 
tremely fast tracker ||, providing a coarse track recon- 
struction based on fast digitization of drift chamber sig- 
nals. Only tracks with a measured transverse momen- 
tum larger than 1.5 GeV/c are considered further. Two 
such tracks, matched to distinct hits in the muon sys- 
tems, are required to pass the Level 1 dimuon trigger. No 
additional requirements are made at Level 2. At Level 
3, a detailed reconstruction is performed and opposite 
sign dimuon events with an invariant mass in the range 
2.7—4.0 GeV/c 2 are accepted and written to tape. Stored 
events are reconstructed using the full set of calibrations. 

The following b hadron decay modes are reconstructed: 
B+ -> J/^jK+, B° -> J/xJjK*°, B° s -> J/i/xf), B a -> 



J/ipK s and A£ — > J/-0A . The daughter particles 
are reconstructed in the decay modes K*° — > K + it~, 
4> — > K + K~, Kg — > 7r + 7r~ and A — > pn~ . Charge 
conjugate modes are included implicitly. To reconstruct 
a given final state we try all possible combinations of 
particle hypotheses, since hadronic particle identification 
capabilities are limited. For a given particle hypothesis 
tracks are corrected for energy loss with the correspond- 
ing mass assigned to the track 0. The correction pro- 
cedure 1 1(1 makes use of the material information in a 
GEANT [ll| description of the CDF detector. Material 
is integrated only at radii larger than that of the recon- 
structed decay vertex of long-lived Kg and A particles. 
High track quality is ensured by requiring at least 20 axial 
and 16 stereo hits in the COT. To ensure a precise mea- 
surement of the b hadron decay vertex, only tracks with 
at least 3 axial SVX hits are considered. The SVX hit re- 
quirement is not applied to the daughter tracks from Kg 
and A which tend to decay outside the silicon tracker. 
A muon is reconstructed from tracks matched to track 
stubs in the muon chambers. 

The mass reconstruction begins by constraining the 
two triggered, oppositely charged muons to a common 
vertex. Candidates with a resulting dimuon mass within 
80 MeV/c 2 of the world average J/ijj mass ^2| are se_ 
lected. A pt threshold of 400 MeV/c is required on 
all tracks, except A daughters, for which all available 
tracks are used. A uniform threshold of 2 GeV/c is im- 
posed on the momentum transverse to the beam direc- 
tion of K + , K*° and </> candidates. Mass windows of 
80 MeV/c 2 , 10 MeV/c 2 and 40 MeV/c 2 around the world 
average masses |l2j are required to select K*°, cf> and 
Kg, respectively. Combinations with a pir mass between 
1.10 GeV/c 2 and 1.13 GeV/c 2 are selected as A candi- 
dates. The Kg and A flight directions are reconstructed 
as the vector connecting the J/ip and the Kg or A ver- 
tices. We require the Kg and A momentum vectors to 
be within 0.25° and 0.57° of their flight direction, re- 
spectively. A cut on the b hadron transverse momentum 
Pt > 6.5 GeV/c is applied. For b hadrons, ct is in the 
range ~ 400 — 500 /tm, where r is the proper lifetime. 
In order to reduce background, the 2-dimensional decay 
length of the b hadron, L xy , defined as L xy = t^T , is 

required to exceed 100 /mi, where X is the vector be- 
tween the production vertex and the decay vertex of the 
b hadron. 

To calibrate the momentum scale of the CDF tracking 
system, three values must be determined: the energy lost 
by a track when passing through the material in the inner 
detector, the radius of the tracker (for the track curva- 
ture measurement) and the strength of the magnetic field 
(for the track curvature-to-momcntum conversion). The 
effect of the tracker radius is indistinguishable from the 
magnetic field strength in the calibration, so we neglect 
the tracker radius and describe the procedure in terms of 
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FIG. 1: The reconstructed mass of dimuons from J/ip decays, 
as a function of transverse momentum. The three lines rep- 
resent various stages of corrections: the solid circles indicate 
no correction, the triangles add the material tuning and the 
squares show all corrections including the magnetic field scale. 



a magnetic field calibration. We use a sample of over 1 
million, inclusive J/ip — > decays to calibrate the track 
energy loss and magnetic field. An underestimate of the 
material results in undercorrected energy loss and intro- 
duces a dependence of the reconstructed dimuon mass on 
Pt- The pt dependence is a signature for inadequate ma- 
terial assessment, as an incorrect value of the magnetic 
field produces a shift in invariant mass independent of 
the pt of the reconstructed particle. The first calibration 
step tunes the amount of material to remove the momen- 
tum dependence. Next the magnetic field is scaled so 
that the reconstructed J/ip — ► /z/x mass agrees with the 
world average. The effect of the calibration steps on the 
momentum dependence of the dimuon mass is illustrated 
in Figure^ Final state radiation in the decay of the J/ip 
leads to an asymmetry in the otherwise Gaussian distri- 
bution of the measured dimuon invariant mass. We use 
a Monte Carlo simulation to correct the resulting bias in 
bins of J/ip Pt during calibration. 

After reconstruction of candidates the mass is ex- 
tracted using an unbinned log- likelihood fit, with the sig- 
nal distribution modeled as a Gaussian. The shape of 
the background is investigated using an inclusive Monte 
Carlo sample of b hadron decays. A detailed detector 
simulation based on GEANT is used. The B + — > J/ip K + 
sample contains significant contributions of partially re- 
constructed B° — ► J/ipK*° — > /i + pT K + ix~ and misre- 
constructed B + — ► J/ip n + decays, which are modeled in 
the background probability distribution function. Par- 
tially reconstructed B° — ► J/ipK*° — > [i + n~K + n~ de- 
cays populate the left shoulder in Figure Et) . Events 
from B + — > J/ip ir + decays appear on the right side of 
the signal peak. The misreconstruction of K*° — > K + ir~ 
due to swapped track assignment of K and 7r is taken 



into account for B° — > J/ip K*° decays. No significant 
contributions are found for the other decay modes. Com- 
parisons between data and fits are shown in Figure El 

The systematic uncertainties are summarized in Ta- 
bles [I] and m The largest systematic uncertainties orig- 
inate from the momentum scale and tracker alignment. 
Deviations from the well-measured world averages in the 
tp' — ► /i + /J~, ip' — ► /i + /i~7r + 7r~ and T — * /i + /J - high 
statistics samples are used to determine the uncertainty 
of the momentum scale. The observed deviations, scaled 
by the Q-value of the respective decay, provide an esti- 
mate of the systematic uncertainty. The second uncer- 
tainty of importance originates from the relative align- 
ment of SVX and COT. It is evaluated by comparing 
mass measurements using the combined tracker informa- 
tion to those using the COT information alone. False 
curvature arises from tracker misalignments and results 
in a curvature offset. The net observed effect is an in- 
crease in momentum for negatively charged tracks and 
a decrease for positively charged tracks. False curvature 
effects cancel, to first order, in charge symmetric sam- 
ples. We derive a parametric correction to remove the 
charge dependence and use the mass shift due to this 
correction as a measure of the systematic error. Uncer- 
tainties due to the vertex fit are evaluated using differ- 
ent mass and pointing constraints in the fit. A corre- 
lation between fluctuations in measured curvature and 
reconstructed vertex position is the source of the resolu- 
tion bias uncertainty. Background systematics are deter- 
mined by varying the background description. In the B + 
case, the mass shift due to inclusion and exclusion of the 
misreconstructed B + — > J/ipir + is assigned as systematic 
uncertainty. In the B case, the systematic uncertainty 
is derived by varying the amount of the reflection contri- 
bution within la of expectation. 

We obtain the following results: 



m(B+) 
m(B°) 
m(B° s ) 
m(Ag) 



5279.10 ± 0.41 (stai) 
5279.63 ± 0.53 (stot) 
5366.01 ± 0.73 (stot) 



5619.7 ± 1.2 



(stat) 



0.36 



(ays) 



MeV/c 2 



0.33 (sys) MeV/c 2 , 
0.33 (sys) MeV/c 2 , 
1.2 (sys) MeV/c 2 . 



These results are in agreement with the current world av- 
erages: m(B+) = 5279.0±0.5MeV/c 2 , m(B°) = 5279.4± 
0.5 MeV/c 2 and m{B° s ) = 5369.6 ± 2.4 MeV/c 2 0. 
Our new A° mass measurement agrees with two of the 
three previous measurements and is in excellent agree- 
ment with CDF's Run I measurement 

[3E1E3 The 

achieved precision is better than the current world aver- 
age of m(Ag) = 5624 ± 9 MeV/c 2 0. 

For the mass differences, most systematic uncertainties 
cancel. The momentum scale uncertainty is scaled down 
to the size of the mass difference. The remaining system- 
atic uncertainty originates from differences in the fit mod- 
els. We minimize systematic effects in m(A°) — m(B°) 
by using the B° — > J/ipKg decay mode, which is topo- 
logically similar to A" — » J/ip A , and where we measure 
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FIG. 2: The invariant mass distribution for J/ip K + , J/rp K°* , J/ip <p and J/tp A candidates. The results of the log-likelihood 
fits are superimposed. The fit probability obtained from a \ 2 test is shown. 



TABLE I: Summary of systematic uncertainties for the B 
meson mass measurements in MeV/c 2 . 



Source 

Tracking 
Momentum scale 
Alignment 
False Curvature 
Vertex Fitting 
Resolution bias 



B ^>J/tpK*° B ± ^JH>K ± B°^J/ip<l> 



0.20 

0.18 

0.02 b 

0.10 

0.13 



0.22 

0.18" 

0.02 

0.10 a 

0.13 



0.20 

0.18° 

0.02 6 

0.10 a 

0.13 



Bkg Systematics 
K-ty swap in K*° 
J/tpn contamin. 



0.06 



0.13 



Total Uncertainty 



0.33 



0.36 



0.33 



"from B° 
'from S± 



m(B°) = 5280.46±0.63 f s t a t)> m agreement with the more 
precise mass determination from the B° — > J/ipK*° de- 
cay mode. The uncertainties are summarized in Table lTTTl 
We obtained the following results for the mass differences: 



ro^B*) - m(B°) = 
m(B°) - m(B°) = 
m(A°) - m(B°) = 339.2 



-0.53 ± 0.67 ± 0.14 MeV/c 2 , 
0.06 MeV/c 2 , 
0.1 MeV/c 2 . 



.38 



± 0.90 ± 
± 1.4 ± 



These are the most precise measurements of m(B®) — 
m(B°) and m(Ag) - m(B°) to date. 

We thank the Fermilab staff and the technical staffs 



TABLE II: Summary of systematic uncertainties for the A° 
mass measurement in MeV/c 2 . The high statistics B° values 
have been used for the A° systematics. 



Source 


B° - J/i>K% 


A° - J/^A° 


Tracking 






Momentum scale 


0.2 


0.2 


Alignment 


1.0 


1.0" 


Vertex Fitting 


0.7 


0.7" 


Total uncertainty 


1.2 


1.2 



"from B 







TABLE III: Summary of systematic uncertainties for the b 
hadron mass differences in MeV/c 2 . 



mass difference 



mom. scale fit model total uncert. 





- 777(5°) 


0.00 


0.14 


0.14 


m{B°) 


- 777(5°) 


0.01 


0.06 


0.06 


m(B°) 


- 777(5*) 


0.01 


0.13 


0.13 


m(A°) 


- 777(5°) 


0.05 




0.05 
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